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ABSTRACT 

The nonlinear instability of the boundary layer on a heated flat plate placed in an on- 
coming flow is investigated. Such flows are unstable to stationary vortex instabilities and 
inviscid traveling wave disturbances governed by the Taylor-Goldstein equation. For small 
temperature differences the Taylor-Goldstein equation reduces to Rayleigh s equation. When 
the temperature difference between the wall and free stream is small the preferred mode of 
instability is a streamwise vortex. It is shown in this case that the vortex, assumed to be of 
small wavelength, restructures the underlying mean flow to produce a profile which can be 
massively unstable to inviscid traveling waves. The mean state is shown to be destabilized 
or stabilized to inviscid waves depending on whether the Prandtl number is less or greater 
than unity. 


iThis research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Our concern is with the nonlinear instability of forced convection boundary layers over 
horizontal heated flat plates and the subsequent secondary instabilities of these flows. In 
a previous paper, Hall and Morris (1992), hereafter referred to as HM, considered the lin- 
ear aspects of longitudinal vortex instabilities. In general it was shown that this convec- 
tive instability develops in a nonparallel manner and cannot be adequately described by a 
quasi-parallel stability theory of the type discussed by for example Wu and Cheng (1976), 
Moutsoglou, Chen and Cheng (1984), Chen and Cheng (1981). At higher values of the con- 
trolling stability parameter, the Grashof number, nonparallel effects can be taken care of in 
a self-consistent asymptotic manner provided that the vortex wavelength is not too large. 
In the present paper we will be concerned with extending the work of HM into the strongly 
nonlinear regime, in addition we will be concerned with the subsequent three-dimensional 
unsteady breakdown of the flow induced by large amplitude vortex structures. 

In addition to longitudinal vortex structures both Tollmien-Schlichting and inviscid waves 
are possible causes of instability in a heated boundary layer. If buoyancy forces are not too 
large, then the inviscid modes are found to satisfy Rayleigh’s equation, otherwise they satisfy 
the Taylor-Goldstein equation. We show that, when buoyancy forces are sufficiently large 
enough to alter the zeroth order inviscid instability problem, they also enter the equations 
for the basic state and cause the mean velocity and temperature fields to be coupled. 

In order to make progress with an analytical solution of the strongly nonlinear instability 
problem associated with longitudinal vortex instabilities, we use the asymptotic structure 
given by Hall and Lakin (1988) in the context of small wavelength Gortler vortices. In this 
limit, which corresponds to the far downstream behavior of a fixed wavelength vortex, the 
vortices are confined to a finite part of the boundary layer and indeed where they exist the 
mean state adjusts to allow them to remain neutral. We show that in forced convectivion 
boundary layers the vortices occupy a region adjacent to the wall and that the thickness of 
this region increases linearly with Grashof number. The results we present are for similarity 
solutions of the nonlinear vortex-mean flow interaction problem but previous experience with 
the related Gortler problem, Hall (1988), Hall and Lakin (1988) suggests that the structure 
of the nonsimilar solutions is not significantly different. 

A significant result which we find is that for Prandtl numbers less than unity favorable 
pressure gradient flows can have inflectional streamwise mean velocity profiles induced by 
vortices. This suggests that a major role of vortex instabilities might be to modify favorable 
pressure gradient flows to make them unstable to the relatively much more dangerous inviscid 
instabilities. We investigate the Rayleigh instability problem for such flows and show that 
vortices induce a class of unstable Rayleigh waves over a large range of wavenumbers. 
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The procedure adopted in the rest of this paper is as follows: in §2 we formulate the 
nondimensional equations governing vortex instabilities in heated boundary layers. In §3 
these equations are solved for small vortex wavelengths and in §4 a similarity solution is 
obtained. Some large Grashof number properties of the similarity solution are discussed in 
§5 where we show the importance of the Prandtl number in determining the shape of the 
induced mean flow. In §6 we discuss the inviscid instability of the vortex states discussed in 
§3, 4, 5 and finally in §7 we draw some conclusions. 

2. Formulation 

Consider the viscous flow over a semi-infinite flat plate with a typical lengthscale L in 
the flow direction. Suppose that the fluid speed at infinity is U 0 u e (x*/L ) where x* measures 
distance along the wall. We define the parameters R and G by 

R = UqLv~\ G = L 3 g(T 0 - T^fiv^R- 312 . (2.1a, b) 

Here v is the kinematic viscosity, g the acceleration due to gravity and ft is the coefficient of 
expansion of the fluid. The temperatures To and Tx, respectively denote a constant reference 
temperature and the fluid temperature in the freestream. We concern ourselves with the limit 
R —* oo with G held fixed and determine the strongly nonlinear vortex flows which occur in 
this limit. 

Suppose next that the wall temperature is given by 

T = T 00 + (T^ - Tq)T{x* j L) (2.2) 

and that we define dimensionless variables (x,y,z) by 

(x,y,z) = {x* L~ l ,y* L~ l R l/2 ,z* L~ l R l/2 ), (2.3) 

where y m ,z* denote distance in the normal and spanwise directions respectively. We define 
a dimensionless velocity field ( u,v,w ) associated with (x,y,z) by writing 

(«*, v*, W) = U 0 (u, R~ l l\, R~ 1/2 w). (2.4) 

The pressure p* is then written in the form 

p‘ = pUl {?„(*) + Mx)/R ' /2 + §?(*,!/) + + . . . j , (2.5) 

where p is a typical fluid density. The fluid density is given by 

p' = p[ 1 + -fi(T* - T^)] (2.6) 
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and we shall make the Boussinesq approximation throughout this work so that changes in 
p can be ignored unless they are multiplied by gravity. The temperature field is then made 
dimensionless by writing 

T* = ^ + (To - TooJT. (2-7) 

In the Boussinesq approximation the continuity equation, momentum equations and energy 
equations are: 


V • u = 0, 


- 1/2 _n—-~ 

R R 


(m • V)u 


(x)R- 1/2 -p x 

(- P y G + GT) - Py 


+ V 2 u, 


(2.8a, b, c) 


-Pz 


(u-V)T = -V 2 T, 
a 

where a is the Prandtl number, V = {d x ,d y ,d 2 ) and V 2 = (d 2 + d 2 z +R~ x d 2 x ). The boundary 
conditions appropriate to (2.8) are 


u = 0, T = T(x), y = 0, 
u -» tte, u> -*• 0, T — > 0, y -» oo. 


(2.9a, b) 


Note here that if the instability is caused by for example nonuniform wall heating the con- 
ditions at the wall should be modified as in HM. In the absence of any instabilities we write 
u = (u,u,0), T = T and the basic state is determined in the limit R — » oo by 


Wj + Vy 
UU X + VUy = 

Py = 

uT x + vT y = 

u = T7 = 0, T = 
u — > u e , r — > 


o, 

tigWex “h ^yy ^Px 

T, 

'£yy 

a 

= o, 

0, y -> oo. 


(2.10a -f) 


Here we have set p 0x = -u e u ex and introduced the buoyancy parameter S = G/R which we 
take to be 0(1). The longitudinal vortex structure of HM first occurs for G = 0 (R°) but our 
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ultimate concern is with the interaction of vortex structures and Taylor-Goldstein inviscid 
waves which exist for G/R ~ 0(R°). Hence it is convenient at this stage to allow for a 
‘buoyancy’ coupling between the x,y momentum equations, even though most of the results 
given in this paper will be for the case 5 = 0. Note also that (2.10c) can be integrated and 
the result substituted into (2.10b) to give 


Uli, 


+ VUy = U e U ex + U 


vv 


+ S f y T x dy. 

Joe 


(2.11) 


This coupling has a significant effect on similarity solutions for the basic state. In particular 
if we choose u e = x n then the wall temperature T must be proportional to x 5n / 2 -i /2 f or a 
similarity solution to exist. 

In HM the linearized version of (2.8) was solved for the case n = 1; it was shown that 
at 0(1) values of G the longitudinal vortex instability induced by heating develops in a 
nonparallel manner. At higher values of G interest centres on the right-hand branch of the 
neutral curve where vortices of wavenumber a are neutral for G ~ a 4 . In fact this high 
wavenumber regime is important because a disturbance of fixed physical wavelength will 
correspond to large values of the local Grashof and wavenumbers at sufficiently large values 
of x. Thus the regime is particularly relevant to experiments where the vortex instability is 
caused by very small background disturbances which must evolve over a significant distance 
before they develop in a nonlinear manner. It is with the latter nonlinear problem which we 
will concern ourselves in this paper. 


3. The Strongly Nonlinear Evolution Equations for Small Wavelength 
Streamwise Vortices 

The nonlinear structure which we develop in this section is based on the related analysis 
of Hall and Lakin (1988) for Gortler vortex flows. The major difference between the problems 
however is that in the convection problem small wavelength neutral vortices occur near to the 
boundary whilst in the Gortler problem the vortices move away from the wall to the position 
in the flow where Rayleigh’s criterion is most violated. This difference leads to significant 
differences in the corresponding nonlinear structures. In fact the heated boundary layer 
has some similarities with the strongly nonlinear Taylor vortex problem, Denier (1992). 
However, the fact that the Taylor instability occurs in a bounded region causes even the 
Taylor and convection problems to differ significantly at high values of the appropriate 
control parameters. 

In the convection problem, small wavelength vortices feel the local strength of the desta- 
bilizing bouyancy force through the vertical gradient of the basic temperature field. In 
general, this gradient is a maximum at the wall so the instability is initiated there. In fact 
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it was shown in HM that in the neutral case the dominant balance is between spanwise 
diffusion and buoyancy effects in the y momentum equation and between spanwise diffusion 
and the convective operator in the energy equation. It is this balance which leads us to the 
appropriate generalization of Hall and Lakin (1988) to the present problem. 

Suppose then that a vortex of wavelength a exists in a region of depth 0(1) adjacent 
to y = 0. Since we have assumed that G/R is 0(1) then we must consider the case when 
a = 0 (R 1 / 4 ). Thus we write 

G = Ga 4 + • • • 


(3.1a, b) 


with G/R = S + • • • . 


In order to recover the situation when buoyancy is not important we must then consider the 


second limit S -> 0. The neutral right hand branch modes for a > 1 have temperature, 
pressure and velocity perturbations 6,p,ii,v,w such that 

-r ~ 0(a 2 ), t ~ 0(1), t ~ 0(a 2 ), -r ~ 0(a). 

Q p u w 


(3.2) 


The first of the above scalings follows from our earlier remarks concerning the balances 
GO ~ v zz , 9 ZZ ~ vT y whilst the remaining balances follow from the choice of balances 
Pz ~ w zz , u zz ~ vu y , v y ~ w z in the z,x momentum equations and the continuity equation 
respectively. We then fix the overall size of the disturbance so that the mean flow and 
temperature field corrections driven by the disturbance are comparable with the unperturbed 
flow. Thus we require that u ■ Vu ~ 0(1) ~ u ■ V$; we therefore expand u,v,w,p and F 
appearing in (2.8) in the form 


u = u Q + a 2 / 3 ui(x, y) H + a l [U 0 (x,y)E + C.C.] H 


v = v 0 + a 2 / 3 vi(x,y) + ■ ■ • + a[Vo(x,y)E + C.C.] -\ 

w = [W 0 {x,y)E + C.C.} + "- (3.3a, b, c, d, e) 


p = [a 2 P 0 (x,y)E + C.C.] d 

T = Tq(x, y) + a~ 2 / 3 Ti(x,y) + • • • + a~ l [0oE + C.C.] + • • • . 

In the above expressions the function 2 E = e iaz , C.C. denotes complex conjugate and ••• 
denotes terms smaller (in terms of a) than those immediately before this symbol. Notice 
also that the mean (in terms of z) part of the pressure must be expanded as 

p = %{ x > y) + a ~ 2 / 3 qi( x , y) h — ( 3 - 3f ) 

Furthermore, we have anticipated that the correction terms in (3.3) to the mean and funda- 
mental are 0(a -2 / 3 ) smaller than the dominant terms. This choice can be inferred from the 
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fact that the correction terms to the mean must be comparable to the depth of the transition 
layer in which the vortices decay to zero. Since the depth of this layer is governed by the 
scalings of the linear neutral problem it must be of depth 0(a~ 2 / 3 ), see Hall and Lakin (1988). 

Equations for the vortex 


If we substitute (3.3) into (2.8) and retain the leading 


order fundamental terms we obtain: 


Voy + iWo — 0 , 

Uq + Vou 0j( = 0, 
G0 o -V o = 0, 
Wo + iPo = 0, 


(3.4a, b, c, d, e) 


— + VoToy = 0. 
a 

These equations determine only U 0 , V 0 , W 0 , and P 0 in terms of 6 0 , and in fact the equations 
are only consistent if 

crGToy + 1 = 0 . 

We integrate this equation to give 


r 0 = T(x) - i, 

CrO* 


(3.4f) 


where we have satisfied the required boundary condition on T 0 at y — 0. The temperature 
profile (3.4f) is that temperature distribution which enables a vortex of wavenumber a > 1 
to remain neutral. Without any loss of generality we will now assume that the fundamental 
temperature disturbance 9 0 is real. 


Equations for the mean : 

If we substitute the expansions (3.3) into (2.8) and equate the leading order mean terms 
in the continuity, x, y momentum equations and the energy equations and use (3.3) we obtain 

u Ox + Voy = 0, 


+ Vo u Qy 


%y 


UqTq x + VqT by 



(3.5a, b, c, d) 
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The y-momentum equation above can be integrated directly to give 

% — f Tody + %(x, 0). (3.5e) 

Jo 


The function T 0 appearing in (3.5c) can then be replaced by 1 - ^ and then (3.5a, b,d) can 
be integrated to find V 2 ,u 0 ,v 0 . It is instructive at this stage to seek a similarity solution of 
(3.5) in order to gain some insight into the nonlinear structure we have obtained. 

We recall that for the Falkner-Skan profile with u e ~ x n buoyancy forces can only be 
retained within the similarity solution structure if T ~ x^. However (3.4c) implies that a 
similarity solution in the presence of vortices is possible only if T ~ z (1_n)/2 - Thus a similarity 
solution in the presence of both buoyancy effects and vortices is possible only if n = 1/3. In 
fact the case n = 1/3 is of particular interest in vortex-wave interaction theory, since, as will 
be remarked upon later, the only possible similarity solution when small wavelength vortices 
interact with Rayleigh or Taylor-Goldstein waves has n = 1/3. In the next section we will 
discuss the n = 1/3 case in detail, before doing so we will complete the description of the 
flowfield for the more general non self-similar case. However we stress that the results which 
we obtain for this particular choice of n are typical of the other similarity solutions. 

The equations (3.5a, b,d) must in general be integrated numerically with % given by 
(3.5c) and T 0 given by (3.4f). However, we can think of (3.5d) as a first order equation in y 
for Vq and we formally write the solution in the form 


R 0 2 = V o 2 (x,0)-2 6 



dy. 


(3.6) 


The function V o (x,0) remains unknown at this stage and for small values of y the integrand 
behaves like u 0y (x,y)yT x so that if we assume there is no reversed flow then this quantity 
is positive if T x > 0. Thus for similarity solutions with n < 1, V<f is a positive decreasing 
function near the wall and we expect that, since we do not anticipate the presence of vortices 
everywhere, the vortex will vanish at some value of the similarity variable T). In fact for 
similarity solutions with values of n not satisfying this inequality the integrand becomes 
positive for large enough T] , thus the only major change is that the maximum of V& occurs 
away from the wall. We also expect that (3.6) will determine V 0 2 in a finite range of valuesfor 
y in the non self-similar case. 

The solution (3.6) is therefore valid for 0 < y < y and more precisely near y, Vo ~ |j/ — y\ 
and viscous effects come back into play. In fact a layer of depth <T 2 ' 3 is required at y in 
order to allow V 0 to adjust from algebraic decay for y - y > a ~ 2 / 3 to exponential decay for 
y -y a -2/3 . The required structure in this layer is virtually identical to that given by 
Hall and Lakin (1988) and so will not be repeated here. It suffices to say that Vo satisfies 
the second Painleve equation in this layer and that across the layer Uo,Uo y ,Vo,To,T 0 y and q Q 
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are continuous. Thus above the transition layer we retain the expansions (3.3) but with all 
the z dependent terms set equal to zero so that the leading order problem for uq,T>o, etc. in 
this upper layer is identical to that for the unperturbed flow. We can therefore write down 
the following ‘composite’ problem for the whole flow field. 

+ ^0 y ~ Oj 

UoU 0x + V 0 u 0y = -Sq Qx (x, 0 ) - S J T 0x dy + u e u ex + u 0yy 

+ \h(v^{u 0 X}„ 

u 0 T 0t + v 0 T 0y = + | HWWoX),, 

H(vS) {r„ - r(x) + 4-} = o, 

where H is the Heaviside function 

f 1, 3 > 0 

H( S ) = . 

.0, a < 0. 

The appropriate boundary conditions are 

To = Vb(x,0),u o = v 0 = 0, T 0 = T, y = 0 ' 

u 0 — ► u e , T 0 — * 0, y — > oo . ( 3 . 8 ) 

and % x (x,0) = - f™T 0x dy. 

We note that at the point y where V 0 2 = 0 the functions uo,u 0 y,Vo,To,T 0 y are to be made 
continuous and, since (u 0 ,T 0 ) — > (u e ,0) as y oo, then Vo(x,0), u 0y (x,0), must adjust 
themselves in order that this limit is achieved. Note also that the position y where Vo 
vanishes will also be a function of x. In fact for numerical reasons it is more convenient to 
treat <?o*( x !0) as a third unknown to be iterated upon on until ? o (x,0) + / 0 °° T 0 dy = 0. The 
numerical free boundary value problem specified by (3.7, 3.8) can be solved in principle by 
adapting the procedure described in Hall and Lakin (1988), however it was found in that 
paper that the form of the non-self-similar solution can be inferred from the self-similar ones 
by varying the parameter corresponding to G. 

It remains for us to discuss the nature of the flow described above in the neighborhood 
of the wall. We recall that the solution which we have obtained has u = v = 0, T Q = T at 
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the wall but V o (^,0) ^ 0. Of course the total flow must satisfy the no-slip conditions at the 
wall so an inner boundary layer must be present as the wall is approached. We can see from 
(3.3a) that in the limit y -> 0 the solution calculated above is such that 

u ~ yu Oy {x,0) + ■ • • + a 1 [Uo{x, Q)E + •■•]' 

which suggests that a new structure will emerge when y falls to 0(a _1 ). Thus we define a 
new variable ( by 

C = a V 

and seek a solution for ( = 0(1). An examination of the higher harmonics in (3.3), which 
are smaller than the fundamental for y = 0(1), shows that for C = 0(1) all modes are 
comparable. This is because the cascade of energy from the fundamental down to the 
harmonics is enhanced to such an extent by the vertical diffusion of vorticity, now comparable 
with spanwise diffusion, that the energy in the different modes is of a similar size. 

Thus in the neighborhood of the wall we must replace (3.3) by 

u = a _1 u(x, C, -?) + 0(a _5 ^ 3 ), 

v = a£>(:r, £, z) + 0(a 1/,:J ), 

w = aw(x, C, z) + 0(a -2 / 3 ), (3.9a, b, c, d, e) 

p = a 2 p + 0(a 4 / 3 ), 

T = T + a- 1 f(x,(,z) + 0{a-^ 3 ), 
whilst p now expands as 

P = q{ x , £) + 0(a -2 / 3 ). 

The zeroth order approximation to (2.8) in the wall layer then becomes, 


+ w z = 0 , 


VU{ + wu z 

— 4" u zz , 

i)V( + wv z 

= GT — Pc + + Vzz9 

VW^ + U)W Z 

= -p z + U><<: + Wzz> 

vT ( + wf z 

= ;{*« + *»}• 


?< = o 


(3.10a, b,c,d,e,f) 
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which must be solved subject to 

u = v — w = 0, 

u -* u 0y (x, 0)C + [U o (x,0)E + C.C . ], ( — ► oo, 

P -» 0, C oo, 

[V o (x,0)E + C.C.], C oo, 

w — ► 0, C — *■ oo, 

9 -» 9 o (z,0),C -» oo. 

and it is easy to show that a solution of the above system can be found by integrating from 
C = oo to ( = 0 using an appropriate asymptotic form for £ » 1. Thus the wall boundary 
layer is passive even though it is fully nonlinear. Hence we can obtain the core solution for 
y = 0(1) without reference to the wall layer problem so we do not discuss further the latter 
problem. We shall now discuss self-similar solutions of the strongly nonlinear interaction 
problem (3.7) - (3.8). 

4. A Self-Similar Solution 

We shall now concentrate on the special case 

u e = a: 1 / 3 , T = x 1 / 3 , 

and introduce the similarity variable rj defined by 

y 

V x 1 / 3 ’ 

We then write 

So = * 1/3 /'M, So = 3 ^tjW' - 2/}, T„ = i‘/ 3 S (,). Vo = VM, ,0 = x 2 ' 3 , «(,). 

The free boundary problem (3.7) - (3.8) can then be written in the simplified ordinary 
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differential equation form: 


r+ l{ 2 /"/ -/« + !} = l Sq (r,)-^q' 

V-j{/'j- 2/s') = -^H(V 2 ){g'V 2 )' 

H^Hg - 1 + -S-) = o, 

«' = 9 

/ = /' = o, T) = 0, /' = 1,9 = 9 = 0, q = oo. 

In addition we require that f,f',q,g and g' are continuous at r/ = rj where V = 0. The 
problem can be solved numerically by making guesses for F(0), q/O), /"(0) and integrating 
(4.1a,b,d) to find V,f, and q with g = 1 - At the point where V = 0 we then set the 
terms proportional to V 2 in (4.1a,b) equal to zero and integrate these equations together with 
(4. Id) to determine f,g and q for rj > rj. This integration is carried out using the computed 
values of g, g' and q found as q — > rj. For arbitrary choices of F(0), /"(0), q/O) the 

conditions /'(oo) = 1, g( oo) = 0, 9(00) = 0 will not be satisfied but we can perform a 
Newton iteration procedure on the wall values of V , f" ,g' until the conditions at infinity are 
satisfied. 

Numerical solutions of (4.1) were in the first instance obtained for a range of values of G 
for S = 0, a = .2, 1,5. In Figure 1 we show the computed values of rj , F 2 (0) and /"( 0). We 
> see that solutions can be obtained for G greater than some finite value, in fact this critical 

value corresponds to the right hand branch of the neutral curve. At the larger values of G 
used in the calculations the results suggest that rj ~ G,F 2 (0) ~ G 2 ,/"(0) ~ G -1 . We shall 
comment further on this asymptotic limit in the next section. 

In Figures 2, 3, 4 we show the functions g and V 2 in terms of rj for different values 

of G with a respectively equal to 1, 5, .2. It can be seen in each case that the boundary layer 
thickens as the Grashof number increases. In Figures 2b, 4b we note the discontinuity in /'" 
which occurs at the transition layer position rj = rj. The discontinuity also occurs for the case 
a — 5 shown in Figure 3b but at the values of G used in the calculations the discontinuity 
is not apparent. In each of the calculations the function V 2 decreases monotonically from 
its value at the wall to zero at g = rj. The temperature profiles shown in Figures 2c, 3c, 4c 
illustrate the large G structure mentioned above. Thus as G increases the interval over which 
g is linear in rj itself increases linearly with G. On the other hand as G is decreased towards 
its linear critical value the temperature profile approaches it’s unperturbed form. It is also 
clear from the calculated velocity and temperature fields that an important consequence of 
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the mean flow being driven by the vortices is that the boundary layer thickness is increased 
from it’s unperturbed value. More precisely we note that in the presence of vortices the 
boundary layer thickness is increased by a factor G from it’s unperturbed value. 

A significant difference between the calculations for the different values of a can be seen 
in Figures 2b, 3b, 4b. We refer to the fact that at the smallest value of a, namely a = .2, the 
function /" has a minimum in the region where the vortices exist and a discontinuity in the 
sign of its derivative across the transition layer. This result is significant because it means 
that the mean velocity profiles associated with /' have inflection points at the minima of 
f" and sign changes in the second derivative of the mean downstream velocity component 
at the transition layers. We stress that no such points were found for the cases a — 1,5, 
this suggests that inflection points can only be created by the vortices below some critical 
value of the Prandtl number. This means that in low Prandtl numbers flows, wall heating 
not sufficiently large to induce Taylor-Goldstein modes because 5 = 0, might still lead to 
highly unstable inviscid Rayleigh waves. The secondary instability of the flows induced by 
the vortices will be discussed in the final section of this paper. 

In Figure 5 we show /'" as a function of t] for the case a — .2. We note that the size of 
the jump in f" across the transition layer increases with G. In addition when G increases we 
see that the region between the wall and the inflection point increases but that f" remains 
relatively small until the transition layer is reached. In order to understand this behavior 
it is instructive to consider the limit G — > oo in (4.1) and see how the underlying flow 
structure evolves. This limit will be discussed in the next section, before doing so we report 
on some calculations we have carried out for the case when buoyancy effects are present. In 
Figures 6a,b,c we show the quantities rj x ,f"{ 0) and t/ 2 (0) for the case a = 1, and 5 = 0, 1, 2. 
The velocity and temperature profiles associated with the 5 = 2 calculations are shown in 
Figures 6d,e,f,g, we note here that the most significant difference between the 5 = 0 and 
5 = 2 results is that buoyancy effects cause f to overshoot it’s free-stream value for a range 
of values of rj. 

5. The Limiting G —* oo Flow Structure 

We stress here that, although our asymptotic analysis is carried out on the similarity 
problem discussed in Section 4, the approach can be applied in a similar way to the nonsimilar 
problem in the limit i->oo with G fixed, see §7 for a brief discussion of this situation. 

In order to develop a large G solution of (4.1) it is first convenient to write (4.1a,c) in 
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the form 


(?)' “ ^{( 6 - -.)/' + 2 /}, 

! {(/“ - 1 - 2//")} + £{(<k - ,)/' + 2/} 


(5.1a, b) 


/"' = 


i + f] 


which are valid if V 2 > 0. 

The numerical results discussed in the previous section suggested that for G 1 the 
vortices are distributed over a region of depth 0(G). Hence we seek a solution of (5.1) with 
rj = 0(G). Actually, since G is always multiplied by o when it appears, it is more convenient 
to treat oG as a large parameter and define 


Go 


(5.2a) 


It then follows that the right and left hand sides of (5.1b) will balance if / ~ Go, V ~ Go, 
we therefore write 

/ = Go^(ip) H , 

(5.2b, c) 

V = GoVty) + ■ - ■ , 

and the zeroth order problem obtained from substituting the above expansions into (5.1) is 


w = {a^-i -2^$ ,/ ) + r / ([i-V’]^ f, + 2^)}y- 2 , 


t ) - 


(5.3a, b) 


If we assume that there is no sublayer structure near ip = 0 then (5.3) must be solved subject 
to 

$ = = 0, V 2 = F 2 (0), ip = 0 (5.4) 

and K(0) is a constant to be determined. It can be shown that for rj > Go, where no vortices 
exist, a large (Go) solution of (4.1) can only be developed if /' — ► 1, g — » 0,rj — ► ( Go ) + so 
in addition to (5.4) we require 


«'= 1, V = 0, Ip = 1. 


(5.5) 


The nonlinear differential system specified by (5.3), (5.4), (5.5) can in general only be solved 
numerically and the solution will fix V(0) and ^"(O). However the expansions (5.2b, c) are 
not uniformly valid as ip — ► 1_ because in that limit V ~ (1 — V 7 ) so f^e constant term in 
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the denominator of (5.1b) will be comparable with V 2 when 1 — ip 0 {Ga) 2 . In order to 
find the appropriate expansions in that layer we observe that in the neighborhood of ip = 1 
we can write 


$ = ip + ip -(- ip 

T = + + 

where ip is a constant determined by the solution in 0 < ip < 1 and ip is given by 


(5.6a, b) 


ip = c{ 1 — ip) 1+a 


(5.7) 


where c is a constant. It follows that in the neighborhood of ip = 1 the expansion (5.2b) 
must be modified to give 


f = Ga 


1 + ip - 


<t> 


+ 


iP(<P) 


{Gay (Gay+>° 


+ 


(5.8) 


where <p = (Gcr) 2 (l — ip). The function ip is then found to be given by 

1+a 

(5. 

The expansion (5.8) is valid until (p = 0 where the vortex vanishes. The solution found above 
is then matched onto the solution of 


ip-- 




2<r(l + ip) 


r+\{ 2/r-/' 2 +i} = o 

satisfying f'(oo) = 1. However, this outer layer is passive and so we do not pursue the 
solution further here. 

It follows from the large Ga solution found above that if a < 1 in this limit u" ~ /'" is 
a maximum in the layer near tj = Ga. In fact from above it can be shown that for <p = 0(1) 

r = -C(l + <,)<,(* - 1) (* + 

Thus f" becomes large if cr < 1, this is consistent with our numerical work. Numerical 
solutions of (5.3) suggest that the constant c is always positive so that f" is positive or 
negative near rj = Ga depending on whether or not a is less or greater than 1. Moreover, in 
the case a < 1, attains its largest value near r) = Ga and this value rv (Ga) 2 ^ 1 a i Since 
we can show that the solution for t] > G 0 must have f" 1 < 0 we see that the results shown 
in Figure 4 which showed a jump in the sign of f" at rj are confirmed and indeed a similar 
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result would be found for any value of a < 1. Furthermore since a small rj solution of (4.1) 
shows that /'" is positive near rj = 0, the continuity of /"' in (0,?/) implies the existence of 
an inflection point in the velocity profile for any Prandtl number less than unity. 

Some numerical solutions of (5.3) are shown in Figures 7, 8, 9. Figure 7 shows $"(0) 
and 1/(0) as a function of a. The computed velocity field in the large Go limit similarly 
agrees well with the full solutions. Of particular interest is the dependence of the function 
on ip and a. In Figures 8a, b we show plots of and T'" as functions of ip and o. The 
results shown confirm that inflectional velocity profiles exist only for a < 1. At finite values 
of Go our calculations suggest again that inflectional profiles exist only for a < 1 though an 
exhaustive check of this possibility in the o-G plane was not carried out. In the case a / 1 

the results of Figure 7 suggest that for any o > 1, is always negative whilst for cr < 1 

is positive in an interval (^>*,1) for some ip* < 1. In order to demonstrate the behavior of 

vjr ,// near <7 = 1 we have in Figure 9 plotted results for the cases o = .98, .99, 1, 1.01, 1.02. We 

see that Vt'" tends to a negative constant as ft — ► 1 for o = 1, whilst for a slightly less than 
1 T'" becomes large and positive in a small interval near ip = 1. For a slightly greater than 
1 the function is always negative but increases significantly in magnitude near ip = l. 
We note that the finite value of for <7 = 1 occurs because the next correction term in 
(5.6a) when cr = 1 is proportional to ip 3 . The implications of the above results for secondary 
instability theory will be discussed in §6. 


6. Rayleigh or Taylor-Goldstein Inviscid Breakdown Induced by Streamwise 
Vortices 

Now let us consider the inviscid instability of the finite amplitude vortex structures 
described in the previous section. The approach we take is based on the work of Hall and 
Horseman (1991) on the instability of Gortler vortices. The most important property of 
an inviscid instability in a boundary layer flow is that it operates on the same streamwise 
lengthscale as the boundary layer thickness. In terms of (2.8) this means that for the inviscid 
wave disturbance £ ~ R} /2 . In addition inviscid waves are time-dependent so that terms 
u t , T t must be added to (2.8b, c). Thus we now write 


^ ti ^ 


f u(x,y,z) \ 
v(x,y,z) 
w{x,y,z) 


U(x,y,z ) \ 

V 



R 1/2 V(x,y ,z) 

w 

— 

R^ 2 W{x,y,z) 

V 


p{x,y,z) 


RP(x,y,z) 

\T) 


T{x,y,z) ) 


0(x,y,z) / 


Here the first term corresponds to the combined mean flow-vortex state driven by buoyancy 
effects whilst the second term represents an inviscid travelling wave disturbance of arbitrarily 
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small amplitude A. We again denote the buoyancy parameter S = G/R and then substi- 
tute (6.1) into (2.8) and equate terms of order A. In the limit R -> oo the zeroth order 
approximation to this system is 


iaV + V y + W z = 0, 
m(u — c)U + Vu y + Wuz = —iaP, 


ia(u — c)V = - Py + Se , (6.2a, b, c, d, e) 

ia(u — c)W = —P z , 
ia(u — c)9 + VT y + WT Z = 0. 

Here we have replaced £ by the wavespeed c and the appropriate boundary conditions are 
^ y — 0, oo. With the function u specified by the steady nonlinear vortex problem 
discussed earlier the system (6.2) and the boundary conditions constitute an eigenvalue 
problem 


a = a(c, S ). 

Note that if the basic state is dependent only on x, y then (6.2) can be reduced to 


qrp 1 

(u — c)(Vyy — Q?V ) — UyyV = : V , (6'3) 

which is the so-called Taylor-Goldstein equation. If buoyancy forces are negligible, 5 = 0, 
(6.3) reduces to Rayleigh’s equation. However when u is a function of z no simple general- 
ization of (6.3) is available, but we can eliminate U, V, W and 6 from (6.2) to give 


f p y 

1 + 1 
V 

( p . ) 

1 „ 

PyTy + P Z T Z 1 

V[u-c] 2 J 



[[u-c] 2 [5r y -a 2 (u- c ) 2 ]_ 


If -S' = 0 the above equation reduces to the equation obtained by Hall and Horseman (1991) 
in their discussion of inviscid instabilities induced by Gortler vortices. The solution of (6.4) 
is of course a nontrivial task bearing in mind the fact that u must in general be determined 
numerically. Here we shall concentrate on the case when S = 0 so that we are in effect 
limiting our analysis to the determination of whether vortex instabilities induced by wall 
heating can trigger a rapidly growing Rayleigh instability. Thus we shall now confine our 
attention to the solution of the eigenvalue problem 


+ JP-\ 

[“-c] Vy [u-c] 2 

Py — 0, y = 0, P ► 0, y — + oo. 


(6.5) 
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with u(x, y i z'j determined by the nonlinear vortex equations in the presence of wall heating 
but with 5 = 0. In particular we will investigate the case when u corresponds to the small 
wavelength solution discussed in §3,4. We recall that in that limit vortices are confined to a 
region of depth 0(1) adjacent to the wall with boundary layers of thickness a , a ' at the 
wall and at the edge of the region of vortex activity respectively. Above the a~ 2/3 transition 
layer the flow is determined by the unperturbed boundary layer equations and the mean 
parts of u , u y with respect to z are continuous across the layer. It follows that a solution 
of (6.5) can be sought with no 2 dependence where vortices are absent. In the lower part of 
the boundary layer, where u expands as in (3.3a), we find that P takes the form 

P = P 0 + a~ 1 P 1 + a' 2 P 2 + a~ 3 P 3 cosaz + a~ 3 P 4 + ■■■ (6.6) 

where P 0 ,Pi etc. depend only on x,y. If we expand a,c in the form 

(a,c) = (a 0 ,co) + a -1 (ai,ci) H (6.7) 

then we find that Po, P 3 satisfy 

( ik \ q2P ° . = 0, (6.8) 

\{u 0 -Co} 2 ) y [uo-co] 2 

and 

P 3 = -2(fio - Co) 2 (UoP 0y {u 0 - co]- 3 )^ . (6.9) 

Equation (6.8) is simply the Rayleigh equation for a unidirectional flow u o so that the vortex 
does not have a direct effect on the zeroth order Rayleigh problem. However it does have a 
significant indirect effect because u 0 is determined by the vortices. In the region above the 
a -2 / 3 transition layer (6.8) again gives the correct zeroth order approximation to the inviscid 
stability problem. Across the layer u o, u q v are continuous whilst uo yy is discontinuous. 
In fact this discontinuity in u 0yy is smoothed out within the transition layer by viscous 
effects. An examination of (6.5) in the transition layer shows that as long as c 0 ^ u 0 in 
this layer then Po, Poy ar e continuous across the layer. Similarly the a -1 wall layer is passive 
so that if we are not concerned with neutral waves propagating downstream with the mean 
downstream fluid speed in the transition layer then it is sufficient for us to solve (6.8) with uo 
determined by the vortices in 0 < y < y and by the boundary layer equations for y > y. For 
definiteness we consider the case when u 0 has the similarity form discussed in §4. Thus we 
write u 0 = with / determined by (4.1), it is convenient to rescale o;o,co by writing 

a 0 = ct Q x ' , Cq = cox ' 
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so that the eigenvalue problem a 0 = a 0 (co,G ) becomes 

f P 0v \ _ _ 

\[f ~ Co] 2 / v [f'-ZoV ~ 

Pa n = 0, T] = 0, P 0 -*■ 0, // -> oo. (6.10) 

In fact the similarity solution we have considered is particularly important in vortex-wave in- 
teraction theory. In the latter theory, see Hall and Smith (1991), a small amplitude Rayleigh 
oi Tollmien-Schlichting wave system interacts with itself to drive a strong vortex field which 
itself acts back on the wavefield. In the present context a Rayleigh-vortex interaction occurs 
if A in (6.1) is choosen appropriately. In this case the wave system drives the vortex in the 
critical layer and the fact that the wave system must remain neutral as it moves downstream 

means that similarity solutions of the inviscid equation describing the wave are possible only 
if n = i. 

In the absence of wall heating the basic state has /' determined by the Falkner-Skan 
problems for a pressure gradient proportional to aT 1 / 3 , and since the velocity profile is 
non-inflectional no unstable inviscid eigenvalues exist. We saw in the previous section that 
inflectional streamwise velocity profiles are generated by wall heating whenever the Prandtl 
number a is less than unity. 

In Figures 10, 11, 12, 13 we show results we have obtained by solving (6.10) for the case 
a &nd a range of values of G. In our calculations we have kept qq real and computed 
the corresponding complex value of c 0 . We see in Figure 10 that at each value of G there is 
a band of unstable modes present to the right of a finite value of Oq . Each mode becomes 
neutral and disappears when the wavespeed is equal to the fluid speed at the inflection point 
of the velocity profile. Note here that all our calculations were for cases where such an 
inflection point exists, at sufficiently small values of G the basic mean profile approaches 
its unperturbed values and no inflection points exist. The unstable mode persists for all q 0 
greater than the neutral value but has (& 0 co),- tending to zero as ao — > oo, in this limit Co r 
approaches the fluid speed at the discontinuity in /"' at the transition layer. An analysis 
of that limit shows that the present analysis breaks down when a 0 = 0(loga). In this 
wavenumber regime the Rayleigh wave has a two layer structure (of depth a -2/3 , (log a) -1 ) 
at the transition layer and is effectively zero elsewhere. We do not give details of the behavior 
foi do ~ | log a | because the mode is neutral there and the most unstable growth rates occur 
for a 0 = 0(1). We see that the wavenumber of the most unstable mode increases as the 
Grashof number increases, this is of course due to the thickening of the boundary layer. 
Some eigenfunctions associated with the modes are shown in Figures 12, 13. We see the 
concentration near fj of the eigenfunction of the second mode as its growth rate approaches 
zero at large d 0 - 
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7. Conclusion 


We have shown that wall heating produces large amplitude streamwise vortex structures 
which completely alter the boundary layer in which they develop. We have concentrated our 
attention on self-similar flows which enabled us to solve the mean flow-vortex interaction 
problem by reducing it to a set of nonlinear differential equations. A significant result which 
we found was that the mean-state modified at zeroth order can have inflection points whereas 
the unperturbed state does not. For the special case when the driving pressure gradient is 
proportional to x -1 / 3 we found that inflectional profiles are only generated when the Prandtl 
number is less than unity. 

The importance of the inflectional profiles is that they are highly unstable to Rayleigh 
waves growing on a streamwise lengthscale 0(i2 -1 / 2 ) shorter than that over which the near 
state develops. This means that a boundary layer inviscidly stable in the absence of wall 
heating can be made massively inviscidly unstable by streamwise vortex restructuring of 
the boundary layer. It is of course relevant to question whether the significantly different 
mean flow character corresponding to the cases cr < 1, cr>lisa function of the particular 
similarity flow considered. In order answer this question we note that the more general form 
of the similarity solution given by (4.1) with S = 0 is 

n — 1 

v = yx 2 , 

„ ... . i"' 1 \n — 1 n + 1 ' 

uo = x f {rj),v o = — nf + — f\ » 

Vo = V(r}),T 0 = x^girj), 

S'" + { n -^U" - n s a + »} = V 7 ), 

1 Z (7.1) 

-g" + ^fg' + ^ f's = -\{g'v 2 )H(v 2 ), 

a l l £ 

H(V’)(g -l + £)=0 

f = f = 0,t) = 0,/' = l,g = 0,r) = oo, 

/) f',g,g' continuous at g = rj, where V(rj) = 0. 

In order to see whether inflectional profiles exist at large Go we can repeat the analysis of §5 
by seeking a solution for o. The expansion procedure follows exactly that as of §5 and the 
key result is that the correction term 4> in (5.6a) is independent of n and again proportional 
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to (—ip + 1) <T+1 . This means that /'" becomes large as ip —> l- and changes sign when a 
passes through 1. This result can be used to infer that /'" is positive as ip — ► 1— for a < 1 
and coupled with the fact that /"' must be negative for r\ > rj and for t] < 1 we find that for 
any n the large Gcr limit leads to inflectional profiles for a < 1. This argument suggests that 
the results obtained for the special case n — 1/3 are typical and that inflectional profiles are 
produced whenever a < 1 at sufficiently large values of G 0 . In fact the large Grashof number 
analysis of §5 can be reformulated as a large x asymptotic solution of the full interactive 
problem. This can be done for u e ~ x n ,T ~ x m for any positive m and it is found that 
inflectional profiles are again only created for Prandtl numbers less than unity. This suggests 
that all heated boundary layer flows of fluids with Prandtl number less than unity are caused 
to become inviscidly unstable by a streamwise vortex restructuring of the flow. 
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Figure la. The quantities /"( 0), V 2 (0), rj for the similarity solution with <7 = 1, 5 = 0. 



Figure lb. The quantities /"( 0), V 2 (0), T) for the similarity solution with <7 = 5, 5 = 0. 








Figure 2c. The function g(r]) for S = 0, a = 1, G = 2., 4., 6., 8. 
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(7 = . 2 , G= 17 . 5 , 27 . 5 , 37 . 5 , 47.5 
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Figure 10. The growth rate d 0 c ot - as a function of d 0 for the case 5 = 0, 
G = 17.5,22.88,24.88,28.88,32.88,36.88,40.88,44.88. 
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